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ABSTRACT 

We study the properties of gas inside and around galaxy haloes as a function 
of radius and halo mass, focussing mostly on z — 2, but also showing some results 
for z — 0. For this purpose, we use a suite of large cosmological, hydrodynamical 
simulations from the Overwhelmingly Large Simulations project. The properties of 
cold- and hot-mode gas, which we separate depending on whether the temperature 
has been higher than 10 55 K while it was extragalactic, are clearly distinguishable 
in the outer parts of massive haloes (virial temperatures ^> 10 5 K). The differences 
between cold- and hot-mode gas resemble those between inflowing and outflowing gas. 
The cold-mode gas is mostly confined to clumpy filaments that are approximately in 
pressure equilibrium with the diffuse, hot-mode gas. Besides being colder and denser, 
cold-mode gas typically has a much lower metallicity and is much more likely to be 
infalling. However, the spread in the properties of the gas is large, even for a given mode 
and a fixed radius and halo mass, which makes it impossible to make strong statements 
about individual gas clouds. Metal-line cooling causes a strong cooling flow near the 
central galaxy, which makes it hard to distinguish gas accreted through the cold and 
hot modes in the inner halo. Stronger feedback results in larger outflow velocities 
and pushes hot-mode gas to larger radii. The gas properties evolve as expected from 
virial arguments, which can also account for the dependence of many gas properties on 
halo mass. We argue that cold streams penetrating hot haloes are observable as high- 
column density H I Lyman-a absorption systems in sightlines near massive foreground 
galaxies. 

Key words: galaxies: evolution - galaxies: formation - galaxies: haloes - intergalactic 
medium - cosmology: theory 



1 INTRODUCTION 

The gaseous haloes around galaxies grow by accreting gas 
from their surroundings, the intergalactic medium, which 
is the main reservoir for baryons. The galaxies themselves 
grow by accreting gas from their haloes, from which they 
can form stars. Some of the gas is, however, returned to 
the circumgalactic medium by galactic winds driven by su- 
pernovae (SNe) or active galactic nuclei (AGN) and by dy- 
namical processes such as tidal or ram pressure forces. Such 
interactions between the different gas phases are essential 
for galaxy formation and evolution. 

The physical state of the gas in and around haloes will 
determine how fast galaxies grow. Quantifying and under- 
standing the properties of the gas is therefore vital for theo- 
ries of galaxy formation. It is also crucial for making predic- 
tions and for the interpretation of observations as the phys- 
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ical state of the gas determines how much light is absorbed 
and emitted. 

Theoretical and computational studies of the accretion 
of gas onto galaxies have revealed the existence of two dis- 
tinct modes. In the first mode the inflowing gas experi- 
ences an accretion shock as it collides with the hot, hy- 
drostatic halo near the virial radius. At that point it is 
shock-heated to temperatures similar to the virial value and 
typically remains part of the hot halo for longer than a dy- 
namical time. If it reaches a sufficiently hi gh density, it can 
cool r adiatively and settle i nto a disc (e.g. Rees fe Ostrikerl 
1 19771 : IWhite fc Reesl 1 19781 : iFall fc Efstathioul 1 198(f) This 
mode is referred to as 'hot-mode accretion ' l|Katz et al]|2003l ; 
iKeres et al.ll2005l ). If, on the other hand, the cooling time of 
the gas is short compared to the dynamical time, which is the 
case for haloes of sufficiently low mass, a hot halo is unable 
to form and the accreting gas will not go through a shock 
near the virial radius. The accretion rate t hen depends on 
the infall rate instead of on the cooling rate dWhite fc Frenkl 
ll99ll ; lBirnboim fc Dekell2003l ; lDekel fc BirnboimhuOfj ). Ad- 
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ditionally, simulations have shown that much of the gas 
enters the halo along dense filaments or in clumps, which 
gives rise to short cooling times, even in the presence of a 
hot, hydrostatic halo. This denser gas does not go through 
an accretion shock near the virial radius and will therefore 
remain cold until it ac cretes onto the c e ntral galaxy or is 



hit by an outflow (e.g. IKeres et al.ll2005l : IDekel et al1l2009l : 



Ivan de Voort et alJl2011af). We refer to this mode a s 'cold- 
mode accretion' (|Katz et al.ll2003l ; IKeres et al]|2005h . 

Hot- and cold-mode accretion play very different 
roles in the formation of galaxies and their gaseous 



haloes (IKeres et al. 20051; lOcvirk et al.l 12008 


: 


Keres et al. 


2009allbl; Brooks et alj |200S|; IDekel et all 12009 


Crain et al. 


20ld: Ivan de Voort et al.l l2011allbl: IPowell et al.l l201ll: 


Faucher-Giguere et al. l201lf). It has been shown that cold- 



mode accretion is more important at high redshift, when 
the density of the Universe is higher. Hot-mode accre- 
tion dominates the fuelling of the gaseous h aloes of high- 
mass systems (halo mass > 10 12 Mq ; e.g. lOcvirk et al.l 
120081 : Ivan de Voort et alJ Uoilal ^. The importance of hot- 
mode accretion is much reduced when c onsidering accretion 
onto galaxies (as opposed to haloes) (|Keres et al.l [2009a; 
Ivan de Voort et al.l l2011al ) . At z ^ 1 all galaxies accrete 
more than half of their material in the cold mode, although 
the contribution of hot-mode accretion is not negligible for 
high-mass haloes. Cold-mode accretion provides most of the 
fuel for star f ormation and shapes the co smic star formation 
rate density ( van de Voort et al.ll20lfbh . 

Ivan de Voort et al.l (|2011al lbh investigated the roles of 
feedback mechanisms on the gas accretion. They found that 
while the inclusion of metal-line cooling has no effect on 
the accretion onto haloes, it does increase the accretion rate 
onto galaxies, because it decreases the cooling time of the 
hot halo gas. Feedback from SNe and AGN can reduce the 
accretion rates onto haloes by factors of a few, but accretion 
onto galaxies is suppressed by up to an order of magnitude. 
The inclusion of AGN feedback is particularly important for 
suppressing hot-mode accretion onto galaxies, because it is 
mainly effective in high-mass haloes and because diffuse gas 
is more susceptible to outflows. 

Hot, hydrostatic halo gas is routinely studied using X- 
ray observations of galaxy groups and clusters and has per- 
haps even been de tected around individual g alaxies (e.g. 
ICrain et al]|2010l l?l; [Anderson fc BregmanlbOllh . As of yet, 
there is no direct observational evidence for cold-mode accre- 
tion, even though there are claims of individual detections 
in Hi absorption based on the low metalli city and prox- 
imity to a galaxy of the ab sorption system |Ribaudo et al.l 
l201ll ; iGiavalisco et alj 1201 lh . Cosmological simulations can 
reproduce the ob served Hi column density distribution 
l|Altav et al.ll201ll ). They show that cold-mode accretion is 
responsible for much of the observed high column density H I 
absorption at z ~ 3. In particular, most of the detected Ly- 
man limit and low column density damped Lyman-q absorp- 
tion may arise in cold accretion streams (|Fumagalli et al.l 
l201ld ; lrande Voort et"aill201ld ). 

It has also been claimed that the diffuse Lyman-a emis- 
sion detected around s ome high-redshift galaxies is powered 
by cold accretion (e.g. [Fardal et al.l l200ll; iDiikstra fc Loebl 
l2009l ; lGoerdt et al.ll20ld ; lRosdahl fc BlaizotlbOllI 'l. but both 
simulations and observations indicate that the emission 
is more likely scattered light from central H 11 regions 



(e.g. Furlanetto et al.l 20051; Faucher-Giguere et al. 20 id ; 
ISteidel et al.ll2010l ; lHaves et al.ll201ll ; iRauch et al.ll201ll) . 

The temperature is, however, not the only difference 
between the two accretion modes. In this paper we use 
the suite of cosmological hydrodynamical simulations from 
the Overwhelmin gly Large Simulations project (OWLS; 
ISchave et al]|20ld ) to investigate other physical properties, 
such as the gas density, pressure, entropy, metallicity, radial 
peculiar velocity, and accretion rate of the gas in the two 
modes. We will study the dependence of gas properties on 
radius for haloes of total mass ~ 10 12 Mq and the depen- 
dence on halo mass of the properties of gas just inside the 
virial radius. Besides contrasting the hot and cold accretion 
modes, we will also distinguish between inflowing and out- 
flowing gas. While most of our results will be presented for 
2 = 2, when both hot- and cold-mode accretion are impor- 
tant for haloes of mass ~ 10 12 Mq, we will also present some 
results for z = 0, which are therefore directly relevant for ob- 
servations of gas around the Milky Way. We will make use 
of the different OWLS runs to investigate how the results 
vary with the efficiency of the feedback and the cooling. 

This paper is organized as follows. The simulations are 
described in Section [5] including the model variations, the 
way in which haloes are identified, and our method for dis- 
tinguishing gas accreting in the hot and cold modes. In Sec- 
tions|3]and|3]we study the radial profiles and the dependence 
on halo mass, respectively. In Section [5] we discuss the differ- 
ence in physical properties between inflowing and outflowing 
gas. We assess the effect of metal-line cooling and feedback 
from SNe and AGN on the gas properties in Section [6] In 
Section[7]we study the properties of gas around Milky Way- 
sized galaxies at z = 0. Finally, we discuss and summarize 
our conclusions in Section [8] 



2 SIMULATIONS 

To investigate the gas properties in and around haloes, 
we make use of sim ulations taken from the OWLS project 
ijSchave et al.|[2O10h . which consists of a large number of cos- 
mological simulations, with varying (subgrid) physics. Here, 
we make use of a subset of these simulations. We first sum- 
marize the reference simulation, from which we derive our 
main results. The other simulations are described in Sec- 
tion 12.11 For a full descrip t ion of the simulations, we refer 
the reader to ISchave et al.l (|2010h . Here, we will only sum- 
marize their main properties. 

We use a modifie d version of gadget-3 (last de- 
scribed in Springcl 2005), a smoothed particle hydrodynam- 
ics (SPH) code that uses th e entropy formulation of SPH 
<|Springel fc Hernquisdl2002l l. which conserves both energy 
and entropy where appropriate. 

All the cosmological simulations used in this work as- 
sume a ACDM cosmology with parameters derived from 
the WMAP year 3 data, Q m = 1 - fi A = 0.238, O b = 
0.041 8, h = 0.73, cr 8 = 0.74, n = 0.951 (|Spergel et all 
120071 ) . These values are co nsistent with the WMAP year 7 
data l|Komatsu et al.ll201lf ). The primordial abundances are 



1 The most significant discrepancy is in as, which is 8 per cent, 
or 2.3<r, lower than the value favoured by the WMAP 7-year data. 



© 2011 RAS, MNRAS 000,1111221 



Gas properties in and around haloes 3 



Table 1. Simulation parameters: simulation identifier, comoving box size (Lbox)> number of dark matter particles (N, the number of 
baryonic particles is equal to the number of dark matter particles), mass of dark matter particles (rriDM)i initial mass of gas particles 
f" 1 ™ )> number of haloes with 10 115 Mq < Af^alo < 10 12,5 Mq, and number of haloes with more than 100 dark matter particles. 



simulation 


-^box 


N 


T71DM 


™initial 
m gas 


number of haloes with 


number of resolved haloes 




(/i- x Mpc) 




(M s ) 


(M ) 


10 11 - 5 M < M halo < 10 12 - 5 M Q 


at z = 2 


L100N512 


100 


512 3 


5.56 x 10 8 


1.19 x 10 8 


4407 (z = 2) 


32167 


L050N512 


50 


512 3 


6.95 x 10 7 


1.48 x 10 7 


518 (z = 2); 1033 (z = 0) 


32663 


L025N512 


25 


512 3 


8.68 x 10 6 


1.85 x 10 6 


59 (z = 2) 


25813 



X = 0.752 and Y = 0.248, where X and Y are the mass 
fractions of hydrogen and helium, respectively. 

A cubic volume with periodic boundary conditions is 
defined, within which the mass is distributed over N 3 dark 
matter and as many gas particles. The box size (i.e. the 
length of a side of the simulation volume) of the simulations 
used in this work are 25, 50, and 100 h^Mpc, with N = 512. 
The (initial) particle masses for baryons and dark matter are 

1-5 >< 107 ( 5o^"m pc ) 3 M q and 7.0 x 10 7 ( 50 ^r Mpc ) 3 M 0) 
respectively, and are listed in Table [1] We use the no- 
tation L***N###, where *** indicates the box size in 
comoving Mpc//i and ### the number of particles per 
dimension. The gravitational softening length is initially 
3.9 ( 50 h-t M pc ) comoving kpc, i.e. 1/25 of the mean 
dark matter particle separation, but we imposed a maxi- 
mum of 1 ( 50 Mpc ) ^ _1 kpc proper. We use simulation 
REF.L050N512 for our main results. The L025N512 simu- 
lations are used for images, for comparisons between simula- 
tions with different subgrid physics, and for resolution tests. 
The L100N512 run is only used for the convergence tests 
shown in the appendix. 

The abundances of eleven elements (hydrogen, helium, 
carbon, nitrogen, oxygen, neon, magnesium, silicon, sulphur, 
calcium, and iron) released by massive stars (type II SNe and 
stellar winds) and intermediate mass stars (type la SNe and 
asymptotic giant branc h stars) are followed as described in 
W iersma et al. | d2009bh. We assume the stellar initial mass 
function (IMF) of IChabrierl ll2003h. ranging fr om 0.1 to 
100 Mq. As described in Wiersma et a l. (2009a), radiative 
cooling and heating are computed element-by-element in the 
pre sence of the cosmic micro wave background radiation and 
the lHaardt fc Madaul <|200lh model for the UV/X-ray back- 
ground from galaxies and quasars. The gas is assumed to be 
optically thin and in (photo)ionization equilibrium. 

Star formation is modelle d according to the recipe of 
ISchave fc Dalla Vecchial (2008). The Jeans mass cannot be 
resolved in the cold, interstellar med ium (ISM), which coul d 
lead to artificial fragmentation (e.g. iBate fc Burkertlll99"7l ). 
Therefore, a polytropic equation of state Pot oc pgis is im- 
plemented for densities exceeding 71h = 0.1 cm -3 , where 
Ptot is the total pressure and p gaa the density of the gas. 
This equation of state makes the Jeans mass, as well as the 
ratio of the Jeans length and the SPH smoothing kernel, 
independent of the density. Gas particles whose proper den- 
sity exceeds nu 0.1 cm" 3 while they have temperatures 
T ^ 10 5 K are moved on to this equation of state and can 
be converted into star particles. The star formation rate per 
unit mass depends on the gas pressure an d is set to repro - 
duce the observed Kennicutt-Schmidt law (|Kennicuttf l 998 ) . 



Table 2. Simulation parameters: simulation identifier, cooling 
including metals (Z cool), initial wind velocity (f w i n d), initial 
wind mass loading (rj), AGN feedback included (AGN). Differ- 
ences from the reference model are indicated in bold face. 



simulation 


Z cool 


^wind 

(kms" 1 ) 


V 


AGN 


REF 


yes 


600 


2 


no 


NOSN.NOZCOOL 


no 








no 


NO Z COOL 


no 


600 


2 


no 


WDENS 


yes 


density dependent 


no 


AGN 


yes 


600 


2 


yes 



Feedback f rom star formation is impleme nted using the 
prescription of lDalla Vecchia fc Schavel l|2008l ). About 40 per 
cent of the energy released by type II SNe is injected locally 
in kinetic form. The rest of the energy is assumed to be 
lost radiatively. Each gas particle within the SPH smooth- 
ing kernel of the newly formed star particle has a proba- 
bility of being kicked. For the reference model, the mass 
loading parameter r\ = 2, meaning that, on average, the to- 
tal mass of the particles being kicked is twice the mass of 
the star particle formed. Because the winds sweep up sur- 
rounding material, the effective mass loading can be much 
higher. The initial wind v e locity is 600 kms" 1 for the refer- 
ence model. ISchave et all l|2010h showed that these param- 
eter values yield a peak global star formation rate density 
that agrees with observations. 



2.1 Variations 

To investigate the effect of feedback and metal-line cooling, 
we have performed a suite of simulations in which the sub- 
grid prescriptions are varied. These are listed in Table [2] 

The importance of metal-line cooling can be demon- 
stated by comparing the reference simulation (REF) to a 
simulation in which primordial abundances are assumed 
when calculating the cooling rates (NOZCOOL). We also 
performed a simulation in which both cooling by metals 
and feedback from SNe were omitted (NOSN.NOZCOOL). 
To study the effect of SN feedback, this simulation can be 
compared to (NOZCOOL). 

In massive haloes the pressure of the ISM is too high 
for winds with velocities of 600 kms -1 to bl ow the gas out 
of the galaxy (|Dalla Vecchia fc Schavel 120081 ). To make the 
winds effective at higher halo masses, the velocity can be 
scaled with the local sound speed, while adjusting the mass 
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loading so as to keep the energy injected per unit stellar 
mass constant at ~ 40 per cent ( WDENS). 

Finally, we have included AGN feedback (AGN). Black 
holes grow via mergers and gas accretion and inject 1.5 per 
cent of the rest-mass energy of the accreted gas into the sur- 
rounding matter in t he form of heat. The m odel is based on 
the one intro duced by Springel et al.l i|2005f l and is described 
and tested in lBooth fc Schavd ( 20091 1 . who also demonstrate 
that the simulation reproduces the observed mass density in 
black holes and the observed scaling relations between black 
hole mass and central stellar velocit y dispersion and between 
black hole mass and stellar mass. iMcCarthv et~afl (|20ld ) 
have shown that model AGN reproduces the observed stel- 
lar mass fractions, star formation rates, and stellar age dis- 
tributions in galaxy groups, as well as the thermodynamic 
properties of the intragroup medium. 



2.2 Identifying haloes 

The first step towards finding gravitationally bound struc- 
tures is to identify dark matter haloes. These can be found 
using a Friends-of-Friends (FoF) algorithm. If the separa- 
tion between two dark matter particles is less than 20 per 
cent of the average separation (the linking length 6 = 0.2), 
they are placed in the same group. Baryonic particles are 
linked to a FoF halo if their nearest d ark matter neighb our 
is in that halo. We then use subfind |Dolag et alJl200Sh to 
find the most bound particle of a FoF halo, which serves as 
the halo centre. In this work we use a spherical overdensity 
criterion, considering all the particles in the simulation. We 
compute the virial radius, R vlr , within which the average 
density agrees with the prediction of the top-hat spherical 
collap se model in a ACDM cosmology (|Brvan fc Norman] 
1 19981 ). At z = 2 this corresponds to a density of p = 169(p). 

We include only haloes containing more than 100 dark 
matter particles in our analysis, corresponding to a mini- 



mum dark matter halo mass of Mhaio = 10 



Iff 



and 



10 8 ' 9 M in the 100, 50, and 25 h~ 1 Mpc box, respectively. 
For these limits our ma ss functions agree very well with the 
ISheth fc Tormenl l| 19991 ) fit. Table [Q lists, for each simula- 
tion of the reference model, the number of haloes with mass 
10 11 ' 5 M Q < Mhaio < 10 12 ' 5 M and the number of haloes 
with more than 100 dark matter particles. 



2.3 Hot- and cold- mode gas 

During the simulations the maximum past temperature, 
Tmax, was stored in a separate variable. The variable was up- 
dated for each SPH particle at every time step for which the 
temperature was higher than the previous maximum past 
temperature. The artificial temperature the particles obtain 
when they are on the equation of state (i.e. when they are 
part of the unresolved multiphase ISM) was ignored in this 
process. This may cause us to underestimate the maximum 
past temperature of gas that experienced an accretion shock 
at densities nu > 0.1 cm" 3 . Ignoring such shocks is, how- 
ever, consistent with our aims, as we are interested in the 
maximum temperature reached before the gas entered the 
galaxy. Note, however, that the maximum past temperature 
of some particles may reflect shocks in outflowing rather 
than accreting gas. 



Another reason why T max may underestimate the true 
maximum past temperature, is that in SPH simulations a 
shock is smeared o ut over a few smoothing len gths, leading 
to in-shock cooling (|Hutchines fc ThomaslEoOol ) . If the cool- 
ing time is of the order of, or smaller than, the time step, 
then the max i mum temperature will be underestimated. 
ICreasev et alJ (|201 lh have shown that a particle mass of 
10 6 Mq is sufficient to avoid numerical overcooling of ac- 
cretion shocks onto haloes, like in our high-resolution sim- 
ulations (L025N512). The Appendix shows that our lower- 
resolution simulations give very similar results. 

Even with infinite resolution, the post-shock tempera- 
tures may, however, not be well defined. Because electrons 
and protons have different masses, they will have different 
temperatures in the post-shock gas and it may take some 
time before they equilibrate through collisions or plasma 
effects. We have ignored this complication. Another effect, 
which was also not included in our simulation, is that shocks 
may be preceded by the radiation from the shock, which may 
affect the temperature evo lution. Disregarding these issues, 
van de Voort et"afl (|2011aT ) showed that the distribution of 
Tmax is bimodal and that a cut at T max = 10 5 ' 5 K natu- 
rally divides the gas into cold- and hot-mode accretion and 
that it produces similar resul ts as studies based on adaptive 
mesh refinement simulations (jOcvirk et al.1120081 ) . This T max 
threshold was c hosen because the cool ing function peaks at 
10 5-5.5 K ( eg |wiersma etafl l2009aT ). which results in a 
minimum in the temperature distribution. Additionally, the 
UV background can only heat gas to about 10 5 K, which 
is therefore characteristic for cold-mode accretion. In this 
work we use the same T max = 10 5 ' 5 K threshold to separate 
the cold and hot modes. 



3 PHYSICAL PROPERTIES: DEPENDENCE 
ON RADIUS 

The gas in the Universe is distributed in a cosmic web 
of sheets, filaments, and haloes. The filaments also affect 
the structure of the haloes that reside inside them or at 
their intersections. At high redshift, cold, narrow streams 
penetrate hot ha l oes and feed galaxies e ff iciently (e.g. 
Keres et all 120051; iDekel fc Birnboiml 120061; lAgertz et all 



20091 ; ICeverino et al.l |2010| ; Ivan de Voort et alj|2011al ). The 
middle row of Fig. [T] shows the overdensity in several haloes 



with different masses, ranging from Mh 



= 10 1 



(left 



panel) to 10 Mq (right panel), taken from the high- 
resolution reference simulation (REF_L025N512) at z — 2. 
Each image is four virial radii on a side, so the physical 
scale decreases with decreasing halo mass, as indicated in 
the middle panels. To illustrate the morphologies of gas that 
was accreted in the different modes, we show the density of 
the cold- and hot-mode gas separately in the top and bot- 
tom panels, respectively. The spatial distribution is clearly 
different. Whereas the cold-mode gas shows clear filaments 
and many clumps, the hot-mode gas is much more spheri- 
cally symmetric and smooth, particularly for the higher halo 
masses. The filaments become broader, relative to the size 
of the halo, for lower mass haloes. In high-mass haloes, the 
streams look disturbed and some fragment into small, dense 
clumps, whereas they are broad and smooth in low-mass 
haloes. Cold-mode accretion is clearly possible in haloes that 
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Log 10 p/<p> 



Figure 1. Gas overdensity at z = 2 in and around haloes with, from left to right, Mh a io = 10 , 10 , 10 , 10 , and 10 5 Mq 
taken from the simulation REF_L025N512. All images show projections of the gas density in cubes of 4 virial radii on a side. The proper 
sizes of the images are indicated in the panels of the middle row. In the middle column all gas was included. In the top (bottom) row 
we have only included cold- (hot-)mode gas, i.e. gas with T max < 10 K (T max ^ 10 J ' 5 K). The filaments, streams, and dense clumps 
consist of gas that has never been heated to temperatures greater than 10 5 ' 5 K. 



are massive enough to have well-developed virial shocks if 
the density of the accreting gas is high, which is the case 
when the gas accretes along filaments or in clumps. 

Fig. [2] shows several physical quantities for the gas in a 
cubic 1 comoving Mpc region, which is about four times 
the virial radius, centred on the 10 12 Mq halo from Fig. [T] 
These properties are (from the top- left to the bottom-right): 
gas overdensity, temperature, maximum past temperature, 
pressure, entropy, metallicity, radial peculiar velocity, radial 
mass flux, and finally the "hot fraction" which we define as 
the mass fraction of the gas that was accreted in the hot 
mode (i.e. that has T max 10 5 ' 5 K). The properties are 
mass- weighted and projected along the line of sight. The 
virial radius is 264 h~ 1 comoving kpc and is indicated by 
the white circles. 

In Fig. [3]we show the same quantities as in Fig. [2] as a 
function of radius for the haloes with 10 11 ' 5 Mq < Mh a io < 
10 12 ' 5 M at z = 2 in simulation REF_L050N512. The black 
curves show the median values for all gas, except for the 
last two panels which show the mean values. The red (blue) 
curves show the median or mean values for hot- (cold-)mode 
gas, i.e. gas with maximum past temperatures above (below) 
10 5 - 5 K. The shaded regions show values within the 16th and 
84th percentiles. Hot-mode gas at radii larger than 27? v ir is 



dominated by gas associated with other haloes and/or large- 
scale filaments. 

All the results we show are weighted by mass. In other 
words, we stacked all 518 haloes in the selected mass range 
using R/R v i r as the radial coordinate. The black curves in 
Fig- E (except for the last two panels) then show the val- 
ues of the corresponding property (e.g. the gas overdensity 
in the top-left panel) that divide the total mass in each ra- 
dial bin in half, i.e. half the mass lies above the curve. We 
have done the same analysis for volume-weighted quanti- 
ties by computing, as a function of radius, the values of 
each property that divides the total volume, i.e. the sum of 
Wg a s/p, in half but we do not show the results. The volume 
is completely dominated by hot-mode gas out to twice the 
virial radius, reaching 50 per cent at 3i? v ir- Even though the 
volume-weighted hot fraction is very different, the proper- 
ties of the gas and the differences between the properties of 
hot- and cold-mode gas are similar if we weigh by volume 
rather than mass. 

We find that the median density of cold-mode gas is 
higher, by up to 1 dex, than that of hot-mode gas and that its 
current temperature is lower, by up to 2 dex, at least beyond 
0.2_R v i r . The hot-mode maximum past temperature is on 
average about an order of magnitude higher than the cold- 
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a 
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Figure 2. From the top-left to the bottom-right: gas overdensity, temperature, maximum past temperature, pressure, entropy, metallicity, 
radial peculiar velocity, radial mass flux (in solar mass per year per proper Mpc 2 ), and hot fraction in a cubic 1 h~ 1 comoving Mpc 
region centred on a halo of M^aio ~ 10 12 Mq at z = 2 taken from the REF_L025N512 simulation. The white circles indicate the virial 
radius. 



mode maximum past temperature. The median pressure of 
the hot-mode gas only exceeds that of the cold mode by a 
factor of a few, but for the entropy the difference reaches 
2.5 dex. For R > i? v ir the gas metallicity in the cold mode 
is lower and has a much larger spread (four times larger at 
i?vir) than in the hot mode. Cold- mode gas is flowing in at 
much higher velocities, by up to 150 kms -1 , and dominates 
the accretion rate at all radii. Below, we will discuss these 
gas properties individually and in more detail. 



3.1 Density 

The halo shown in Fig. [2] is being fed by dense, clumpy fil- 
aments as well as by cooling diffuse gas. The filaments are 
overdense for their radius, both inside and outside the halo. 
The top-left panel of Fig. shows that the overdensity of 
both hot- and cold-mode gas increases with decreasing ra- 
dius, from ~ 10 at 10i? v ir to ~ 10 2 at i? v ir and to 10 4 at 
0.1-Rvir- The median density of cold-mode gas is higher by 
up to an order of magnitude than that of hot-mode gas for 
all radii 0.1i? v ir < R < 4i? v i r . The cold-mode gas densi- 
ties exhibit a significant scatter of about 2 dex, as opposed 



to about 0.4 dex for hot-mode gas at i? v ir, which implies 
that the cold-mode gas is much clumpier. Beyond 4i? v i r the 
median hot-mode density becomes higher than for the cold 
mode, because there the hot-mode gas is associated with 
different haloes and/or large-scale filaments, which are also 
responsible for heating the gas. 



3.2 Temperature 

Hot gas, heated either by accretion shocks or by SN feed- 
back, extends far beyond the virial radius (top-middle of 
Fig. [2}. Most of the volume is filled with hot gas. The 
location of cold gas overlaps with that of dense gas, so 
the temperature and density are anti-correlated. This anti- 
correlation is a result of the fact that the cooling time de- 
ceases with the gas density. 

For R > 0.2i? v i r the temperatures of the hot- and cold- 
mode gas do not vary strongly with radius (top middle panel 
of Fig. [5J . Note that this panel shows the current tempera- 
ture and not the maximum past temperature. Gas accreted 
in the hot mode has a temperature ~ 10 6 K at R > 0.2_R v i r , 
which is similar to the virial temperature. The median tem- 
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Figure 3. Properties of gas in haloes with 10 Mq < Mhaio < 10 12 ' 5 Mq at z = 2 as a function of radius for all (black curves), 
hot-mode (red curves), and cold-mode (blue curves) gas. Results are shown for the simulation REF_L050N512. Shaded regions show 
values within the 16th and 84th percentiles, i.e., the ±1<t scatter around the median. From the top-left to the bottom-right, the different 
panels show the mass-weighted median gas overdensity, temperature, maximum past temperature, pressure, entropy, metallicity, radial 
peculiar velocity, the mean accretion rate, and the mean mass fraction of hot-mode gas, respectively. 



perature of the hot-mode gas increases slightly from « 2i? vir 
to « 0.2i? v ir because the hot gas is compressed as it falls in. 
Within 0.5-Rvir the scatter increases and around 0.2/Z v i r the 
median temperature drops sharply to ~ 10 4 K. The dra- 
matic decrease in the temperature of the hot-mode gas is a 
manifestation of the strong cooling flow that results when 
the gas has become sufficiently dense to radiate away its 
thermal energy within a dynamical time. The median tem- 
perature of cold-mode gas peaks at slightly below 10 5 K 
around 2R v \ r and decreases to ~ 10 4 K at 0.1 Rvir- The peak 
in the temperature of the cold-mode gas is determined by 
the interplay between photo-heating by the UV background 
and radiative cooling. The temperature difference between 
the two accretion modes reaches a maximum of about 2 dex 
at 0.3-Rvir and vanishes around 0.1-Rvir- 

3.3 Maximum past temperature 

The maximum past temperature (top-right panel of Fig. 
is by definition at least as high as the current temperature, 
but its spatial distribution correlates well with that of the 
current temperature. As shown by Fig. [3] the difference be- 
tween maximum past temperature and current temperature 
is small at R > 7? v ir, but increases towards smaller radii and 
becomes 1 dex for cold-mode gas and 2 dex for hot-mode 
gas at 0.1-Rvir- 

While the temperature of the cold-mode gas decreases 
with decreasing radius, its maximum past temperature stays 



constant at T max ~ 10 5 K. This value of T max is reached 
around 2R v i r as a result of heating by the UV background. 
Both the current and the maximum past temperature of 
the hot-mode gas increase with decreasing radius for R > 
0.3-Rvir- The fact that T max decreases below 0.37? v ir shows 
that it is, on average, the colder part of the hot-mode gas 
that can reach these inner radii. If it were a random subset of 
all the hot-mode gas, then T ma x would have stayed constant. 

3.4 Pressure 

As required by hydrostatic equilibrium, the gas pressure gen- 
erally increases with decreasing radius (middle-left panels of 
Figs. [2] and [3]). However, the median pressure profile (Fig. [3} 
does show a dip around 0.2 — 0.37? v ir that reflects the sharp 
drop in the temperature profiles. Here catastrophic cooling 
leads to a strong cooling flow and thus a breakdown of hy- 
drostatic equilibrium. 

Comparing the pressure map with those of the den- 
sity and temperature, the most striking difference is that 
the filaments become nearly invisible inside the virial ra- 
dius, whereas they stood out in the density and tempera- 
ture maps. However, beyond the virial radius the filaments 
do have a higher pressure than the diffuse gas. This suggests 
that pressure equilibrium is quickly established after the gas 
accretes onto the haloes. 

Fig. [3] shows that the difference between the pressures 
of the hot- and cold-mode gas increases beyond 2R V1I . This 
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is because at these large radii the hot-mode gas is associ- 
ated with other haloes and/or large-scale filaments, while 
the cold-mode gas is intergalactic, so we do not expect them 
to be in pressure equilibrium. Moving inwards from the 
virial radius, the median pressure difference increases un- 
til it reaches about an order of magnitude at 0.3-R v ir- At 
smaller radii the pressures become nearly the same because 
the hot-mode gas cools down to the same temperature as 
the cold-mode gas. 

Although it takes some time to reach pressure equilib- 
rium if the hot gas is suddenly heated, we expect the hot and 
cold gas to be approximately in equilibrium inside the halo, 
because a phase with a higher pressure will expand, lower- 
ing its pressure, and compressing the phase with the lower 
pressure, until equilibrium is reached. While the pressure 
distributions do overlap, there is still a significant difference 
between the two. This difference decreases somewhat with 
increasing resolution, because the cold gas reaches higher 
densities and thus higher pressures, as is shown in the ap- 
pendix. From the example pressure map (Fig. [5| we can see 
that the filaments inside the halo are in fact approximately 
in pressure equilibrium with the diffuse gas around them. At 
first sight this seems at odds with the fact that the median 
pressure profiles are different. However, the pressure map 
also reveals an asymmetry in the pressure inside the halo, 
with the gas to the left of the centre having a higher pressure 
than the gas to the right of the centre. Because there is also 
more hot-mode gas to the left, this leads to a pressure differ- 
ence between the two modes when averaged over spherical 
shells, even though the two phases are locally in equilibrium. 
The asymmetry arises because the hot-mode gas is a space- 
filling gas and the now has to converge towards the centre 
of the halo, which increases its pressure. The cold-mode gas 
is not space filling and therefore does not need to compress 
as much. 



3.5 Entropy 

We define the entropy as 



S 



P(Mm H ) 5/3 



(1) 



where p is the mean molecular weight, ran is the mass of a 
hydrogen atom, and &b is Boltzmann's constant. Note that 
the entropy remains invariant for adiabatic processes. In the 
central panel of Fig. [2] we clearly see that the filaments have 
much lower entropies than the diffuse gas around them. This 
is expected for cold, dense gas. 

For R > i? v i r the median entropy of hot-mode gas is 
always higher than that of cold- mode gas (Fig. [3]). While 
the entropy of the cold-mode gas decreases smoothly and 
strongly towards the centre of the halo, the entropy of the 
hot-mode gas decreases only slightly down to 0.2f? v i r after 
which it drops steeply. Cold-mode gas cools gradually, but 
hot-mode gas cannot cool until it reaches high enough den- 
sities, which results in a strong cooling flow. 

3.6 Metallicity 

The middle-right panels of Figs. [2] and [3] show that the cold- 
mode streams have much lower metallicities than the diffuse, 
hot-mode gas, at least for R > R v i r . The cold-mode gas also 



has a much larger spread in metallicity, four orders of mag- 
nitude at _R v ir, as opposed to only one order of magnitude 
for hot-mode gas. 

The gas in the filaments tends to have a lower metallic- 
ity, because most of it has never been close to a star-forming 
region, nor has it been aff ected by galactic w inds, which 
tend to avoid the filaments ()Theuns et alJl2002h . The radial 
velocity image in the bottom-left panel of Fig. [2] confirms 
that the winds take the path of least resistance. The cold 
mode also includes dense clumps, which show a wide range 
of metallicities. If the density is high enough for embedded 
star formation to occur, then this can quickly enrich the en- 
tire clump. The enhanced metallicity will increase its cooling 
rate, making it even more likely to accrete in the cold mode 
(recall that our definition of the maximum past temperature 
ignores shocks in the ISM). On the other hand, clumps that 
have not formed stars remain metal-poor. The metallicity 
spread is thus caused by a combination of being shielded 
from winds driven by the central galaxy and exposure to 
internal star formation. 

At J?vir the median metallicity of the gas is subsolar, 
Z ~ 1CT 1 Zq for hot-mode gas and Z ~ 1CT 2 Zq for cold- 
mode gas. However, we caution the reader that the median 
cold-mode metallicity is not converged with numerical reso- 
lution (see the Appendix) and could in fact be much lower. 
The metallicity increases towards the centre of the halo and 
this increase is steeper for cold-mode gas. The metallicity 
difference between the two modes disappears at R ~ 0.5-R v ir, 
but we find this radius to move inwards with increasing res- 
olution (see the Appendix). Close to the centre the hot gas 
cools down and ongoing star formation in the disc enriches 
all the gas. The scatter in the metallicity decreases, espe- 
cially for cold-mode gas, to ~ 0.7 dex. 

As discussed in detail bv lWiersma et all (|2009bT ). there 
is no unique definition of metallicity in SPH. The metallic- 
ity that we assign to each particle is the ratio of the metal 
mass density and the total gas density at the position of the 
particle. These "SPH-smoothed abundances" were also used 
during the simulation for the calculation of the cooling rates. 
Instead of using SPH-smoothed metallicities, we could, how- 
ever, also have chosen to compute the metallicity as the ratio 
of the metal mass and the total gas mass of each particle. Us- 
ing these so-called particle metallicities would sharpen the 
metallicity gradients at the interfaces of different gas phases. 
Indeed, we find that using particle metallicities decreases the 
median metallicity of the metal-poor cold mode. For the hot 
mode the median particle metallicity is also lower than the 
median smoothed metallicity, but it increases with resolu- 
tion, whereas the cold-mode particle metallicities decrease 
with resolution. 

While high-metallicity gas may belong to either mode, 
gas with metallicity < 10 -3 Zq is highly likely to be part 
of a cold flow. This conclusion is strengthened when we in- 
crease the resolution of the simulation or when we use par- 
ticle rather than SPH-smoothed metallicities. Thus, a very 
low metallicity appears to be a robust way of identifying 
cold-mode gas. 



3.7 Radial velocity 

The radial peculiar velocity is calculated with respect to 
the halo centre after subtracting the peculiar velocity of the 
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halo. The peculiar velocity of the halo is calculated by taking 
the mass-weighted average velocity of all the gas particles 
within 10 per cent of virial radius. Note that the Hubble flow 
is not included in the radial velocities shown. It is unimpor- 
tant inside haloes, but is about a factor of two larger than 
the peculiar velocity at 10-R v ir- 

The bottom- left panel of Fig. [2] shows that gas outside 
the haloes is, in general, moving towards the halo (i.e. it has 
a negative radial velocity). Within the virial radius, however, 
more than half of the projected area is covered by outflowing 
gas. These outflows are not only caused by SN feedback. 
In fact, simulations without feedback also show significant 
outflows (see Figs. [7] and [SJ. A comparison with the pressure 
map shows that the outflows occur in the regions where the 
pressure is relatively high for its radius (by 0.1 — 0.4 dex, 
see Figs. [S] and |5J). The inflowing gas is associated with the 
dense, cold streams, but the regions of infall are broader 
than the cold filaments. Some of the hot-mode gas is also 
flowing in along with the cold-mode gas. These fast streams 
can penetrate the halo and feed the central disc. At the same 
time, some of the high temperature gas will expand, causing 
mild outflows in high pressure regions and these outflows are 
strengthened by SN-driven winds. 

The hot-mode gas is falling in more slowly than the 
cold-mode gas or is even outflowing (bottom-left panel of 
Fig. |3]) . This is expected, because the gas converts its ki- 
netic energy into thermal energy when it goes through an 
accretion shock and because a significant fraction of the hot- 
mode gas may have been affected by feedback. For hot-mode 
gas the median radial velocity is closest to zero between 
0.3-Rvir < R < 1-Rvir- Most of it is inflowing at smaller radii, 
where the gas temperature drops dramatically, and also at 
larger radii. 

The cold-mode gas appears to accelerate to —150 km 
s _1 towards i? v ir (i.e. radial velocities becoming more nega- 
tive) and to decelerate to —30 km s _1 from i? v i r towards the 
disc. We stress, however, that the behaviour of individual gas 
elements is likely to differ significantly from the median pro- 
files. Individual, cold gas parcels will likely accelerate until 
they go through an accretion shock or are hit by an out- 
flow, at which point the radial velocity may suddenly vanish 
or change sign. If this is more likely to happen at smaller 
radii, then the median profiles will show a smoothly decel- 
erating inflow. Finally, observe that while there is almost 
no outflowing cold-mode gas around the virial radius, close 
to the central galaxy (R < 0.3i? v ir) a significant fraction is 
outflowing. 

3.8 Accretion rate 

The appropriate definition of the accretion rate in an ex- 
panding Universe depends on the question of interest. Here 
we are interested in the mass growth of haloes in a comov- 
ing frame, where the haloes are defined using a criterion that 
would keep halo masses constant in time if there were no pe- 
culiar velocities. An example of such a halo definition is the 
spherical overdensity criterion, which we use here, because 
the virial radius is in that case defined as the radius within 
which the mean internal density is a fixed multiple of some 
fixed comoving density. 

The net amount of gas mass that is accreted per unit 
time through a spherical shell S with comoving radius x = 



R/a, where a is the expansion factor, is then given by the 
surface integral 

M sas (x) = - f a 3 pi^f (2) 
Js a 

= — / pvv&ddS, (3) 
Js 

where a 3 p and dS/a 2 are a comoving density and a comoving 
area, respectively, and the radial peculiar velocity is n ra d = 
ax. We evaluate this integral as follows, 

M gas (R) = - — 5-! -A shelh (4) 

R^KR+dR shc11 

where 

Vahoii = y((R + dR) 3 -i? 3 ), (5) 

A shcll = 4Tv(R+^dRf, (6) 

n is the radius of particle i, and dR is the bin size. Note 
that a negative accretion rate corresponds to net outflow. 

The mass flux map shown in the bottom-middle panel 
of Fig. [5] is computed per unit area for each pixel as 
Simgas.iWrad.i/Vpix, where Vpi x is the proper volume of the 
pixel. The absolute mass flux is highest in the dense fila- 
ments and in the other galaxies outside R v i r , because they 
contain a lot of mass and have high inflow velocities. 

The gas accretion rate M gas (R/-Rvir) computed using 
equation Q is shown as the black curve in the bottom- 
middle panel of Fig. [3] The accretion rate is averaged 
over all haloes in the mass bin we are considering here 
(10 115 M Q < A/h a io < 10 12 5 M Q ). Similarly, the red and 
blue curves are computed by including only hot- and cold- 
mode particles, respectively. The accretion rate is positive at 
all radii, indicating net accretion for both modes. The inflow 
rate is higher for the cold mode even though the hot-mode 
gas dominates the mass budget around the virial radius (see 
Section [3]9}. The hot-mode gas accretion rate is a combina- 
tion of the density and the radial velocity of the hot-mode 
gas, as well as the amount of mass in the hot mode. Gas 
belonging to the cold mode at R > R v i r may later become 
part of the hot mode after it has reached R < R v i T . 

The extended halo is not in a steady state, because 
the accretion rate varies with radius. Moving inwards from 
10i?vir to -Rvir, the net rate of infall drops by about an or- 
der of magnitude. This implies that the (extended) halo is 
growing: the flux of mass that enters a shell from larger radii 
exceeds the flux of mass that leaves the same shell in the di- 
rection of the halo centre. This sharp drop in the accretion 
rate with decreasing radius is in part due to the fact that 
some of the gas at -R > R v ir is falling towards other haloes 
that trace the same large-scale structure. 

Within the virial radius the rate of infall of all gas and of 
cold-mode gas continues to drop with decreasing radius, but 
the gradient becomes much less steep (din M / din R « 0.4), 
indicating that the cold streams are efficient in transporting 
mass to the central galaxy. For the hot mode the accretion 
rate only flattens at R < 0.4R v ir around the onset of catas- 
trophic cooling. Hence, once the hot-mode gas reaches small 
enough radii, its density becomes sufficiently high for cool- 
ing to become efficient, and the hot-mode accretion becomes 
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efficient too. However, even at 0.1-Rvh- its accretion rate is 
still much lower than that of cold-mode gas. 

3.9 Hot fraction 

Even though the average hot fraction, i.e. the mean frac- 
tion of the gas mass that has a maximum past temperature 
greater than 10 5 ' 5 K, of the halo in the image is close to 
0.5, few of the pixels actually have this value. For most pix- 
els /hot is either close to one or zero (bottom-right panel of 
Fig. [2]), confirming the bimodal nature of the accretion. 

The bottom-right panel of Fig. [3] shows that the hot 
fraction peaks around the virial radius, where it is about 
70 per cent. Although the hot fraction decreases beyond the 
virial radius, it is still 30 per cent around 10i? v ir. The hot- 
mode gas at very large radii is associated with other haloes 
and/or large-scale filaments. Within the halo the hot frac- 
tion decreases from 70 per cent at R v i T to 35 per cent at 
0.1i?vir. While hot-mode accretion dominates the growth of 
haloes, most of the hot-mode gas does not reach the centre. 
Cold-mode accretion thus dominates the growth of galaxies. 



4 DEPENDENCE ON HALO MASS 

In Fig. U we plot the same properties as in Fig. [3] as a func- 
tion of halo mass for gas at radii 0.8i? v ir < R < R v i T , where 
differences between hot- and cold-mode gas are large. Grey, 
dashed lines show analytic estimates and are discussed be- 
low. The dotted, grey line in the top-left panel indicates 
the star formation threshold, i.e. na = 0.1 cm -3 . The dif- 
ferences between the density and temperature of the hot- 
and cold-mode gas increase with the mass of the halo. The 
average temperature, maximum past temperature, pressure, 
entropy, metallicity, absolute radial peculiar velocity, abso- 
lute accretion rate, and the hot fraction all increase with 
halo mass. 

We can compare the gas overdensity at the virial ra- 
dius to the density that we would expect if baryons were 
to trace the dark ma tter, p v ir- We assume an NFW profile 
l|Navarro et al.l 1996). take the mean internal density rela- 
tive to the critical density at redshift z, A c (p), fro m spher- 
ical collapse calculations (I Bryan fc Norman 19981) and the 
halo mass-concentration relation from IPuffv et all (|2008l ) 
and calculate the mean overdensity at i? v ir- This is plot- 
ted as the dashed, grey line in the top-left panel of Fig. [4] It 
varies very weakly with halo mass, because the concentra- 
tion depends on halo mass, but this is invisible on the scale 
of the plot. For all halo masses the median density is indeed 
close to this analytic estimate. While the same is true for the 



hot-mode gas, for high-mass haloes (M ha i > 10 12 M ) the 
median density of cold-mode gas is significantly higher than 
the estimated density and the difference reaches two orders 
of magnitude for Mhaio ~ 10 13 Mg. A significant fraction 
of the cold-mode gas in these most massive haloes is star 
forming and hence part of the ISM of satellite galaxies. The 
fact that cold-mode gas becomes denser and thus clumpier 
with halo mass could have important conseque nces for the 
formation of clumpy ga l axies at high redshift |Dekel et al.l 
120091 ; lAgertz et all 1 20091 ; ICeverino et al.ll2010l ). 

The blue curve in the top-middle panel shows that the 
median temperature of the cold-mode gas at 7? v ir decreases 



slightly with halo mass, from 40,000 K to 15,000 K. This 
reflects the increase in the median density of cold-mode gas 
with halo mass, which results in shorter cooling times. The 
median temperature of the hot-mode gas increases with halo 
mass and is approximately equal to the virial temperature 
for Mhaio > 10 J M©. The virial temperature is plotted as 
the dashed, grey line and is given by 



G 2 //oO m A, 
54 

9.1 x 10 5 K 



)" S =A4^(i + -), 

M halo 



10 12 M Q 



2/3 



1 + z 



(7) 
(8) 



where G is the gravitational constant, Ho the Hubble con- 
stant and n is assumed to be equal to 0.59. 

While much of the gas accreted onto low-mass haloes in 
the hot mode has a temperature smaller than 10 5 ' 5 K and 
has therefore already cooled down substantially, there is 
very little overlap in the current temperatures of gas ac- 
creted in the two modes for haloes with T V - 1T > 10 6 K. 
Because the cooli ng rates decrease wit h temperature for 
T > 10 5 - 5 K (e.g. IWiersma et~alll2009al ). most of the hot- 
mode gas in haloes with higher temperatures stays hot. For 
the same reason, the median temperature of all gas rises 
sharply at M « 10 11 ' 5 M (T vSr « 10 5 5 K) and is roughly 
equal to T vir for M talo > 10 12 M . 

The top-right panel shows that the median maximum 
past temperature of gas at the virial radius is close to the 
virial temperature, which is again shown as the grey dashed 
curve, for the full range of halo masses shown. Some of the 
gas does, however, have a maximum past temperatures that 
differs strongly from the virial temperature. The largest dif- 
ference is found for cold-mode gas in high-mass haloes. Be- 
cause of its high density, its cooling time is short and the 
gas does not shock to the virial temperature. The maximum 
past temperature of gas accreted in the cold mode is close 
to 10 5 K for all halo masses. For Af ha i < 10 10 5 M Q this 
temperature is higher than the virial temperature. The gas 
in low-mass haloes has not been heated to its maximum 
temperature by a virial shock, but by the UV background 
radiation or by shocks from galactic winds. Heating by the 
UV background is the dominant process, because simula- 
tions without supernova feedback show the same result (see 
Fig. [SJ). The maximum past temperature of hot- mode gas 
follows the virial temperature closely for high-mass haloes. 
For T vir < 10 5 5 K (Mhaio < 10 11,5 M ) the maximum past 
temperature of the hot-mode gas remains approximately 
constant, at around 10 5,r K, because of our definition of 
hot-mode gas (T max > 10 5 5 K). 

The pressure of the gas increases roughly as M 2 ^ 3 
(middle-left panel). We can estimate the pressure at the 
virial radius from the virial temperature and the density 
at the virial radius. 



Pvir 



^virpvir 

iim H 



(9) 



where /j, is assumed to be equal to 0.59. This pressure is 
shown by the dashed, grey line. The actual pressure is very 

2 Note that for haloes with T v ; r < 10 5 5 K the median hot-mode 
temperature is affected by the requirement T max > 10 J ' 5 K (our 
definition of hot-mode gas). 
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Figure 4. Properties of gas at 0.8-R v i r < R < R vlI at z = 2 as a function of halo mass for all (black curves), hot-mode (red curves), 
and cold-mode (blue curves) gas. Results are shown for the simulation REF_L050N512. Shaded regions show values within the 16th 
and 84th percentiles, i.e., the ±l<x scatter around the median. From the top-left to the bottom-right, the different panels show the 
mass-weighted median gas overdensity, temperature, maximum past temperature, pressure, entropy, metallicity, radial peculiar velocity, 
the mean accretion rate, and the mean mass fraction of hot-mode gas, respectively. The horizontal, dotted line in the first panel indicates 
the threshold for star formation (njj = 0.1 cm -3 ). The dashed, grey curves show analytic estimates from virial arguments. 



close to this simple estimate. It scales with mass as the virial 
temperature because the density at the virial radius is nearly 
independent of the mass. For all halo masses the median 
pressure of the gas accreted in the hot-mode is about a factor 
of two higher than the median pressure of the cold-mode gas. 

The central panel shows that the entropy difference be- 
tween hot- and cold-mode gas increases with halo mass, be- 
cause the entropy of hot-mode gas increases with halo mass, 
whereas the cold-mode entropy decreases. The hot-mode gas 
follows the slope of the relation expected from virial argu- 
ments, 



5vir 



Pvir(^H) 5/3 
1 5/3 



(10) 



where /j, is assumed to be equal to 0.59. S v h- is shown as the 
dashed, grey line. 

The middle-right panel shows that the median gas 
metallicity at the virial radius increases from ~ 10 -2 Zq 
for M ha io ~ 10 10 M Q to ~ 10 _1 Zq for 10 13 M©. This in- 
crease reflects the increased fraction of hot-mode gas (see the 
bottom-right panel) and an increase in the median metal- 
licity of the cold-mode gas, which is probably due to the 
fact that a greater fraction of the gas resides in the ISM of 
satellite galaxies for more massive haloes (see the top-left 
panel). The scatter in the metallicity of the cold-mode gas 
is always very large. The hot-mode gas has a median metal- 



licity ~ 10" 1 Z Q for all halo masses, which is similar to the 
predicted metallicity of the warm- hot intergalactic medium 
(|Wiersma et al.ll201ll ). 

The black curve in the bottom-left panel shows that for 
all halo masses more mass is falling into the halo than is 
flowing out. The radial velocity distributions are, however, 
very broad. A substantial fraction of the hot-mode gas, more 
than half for Mhaio < 10 11 ' 5 M , is outflowing at R vil . Cold- 
mode gas is predominantly inflowing for all masses, but the 
fraction of outflowing gas becomes significant for Mhaio < 
10 11 ' 5 M Q . 

As expected, the gas at the virial radius falls in faster for 
higher-mass haloes and the absolute velocities are larger for 
cold-mode gas. We can compare the radial peculiar velocity 
to the escape velocity, 



2GM ha 



(ii) 



where we assumed a matter-dominated Universe and used 



-Rvir 



2GM 



1/3 



H%n m A c j i + z 

114 kpc 



M halo V /3 fl + z"- " 
10 12 Mp 



(13) 
(14) 
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We show — « csc by the dashed, grey curve. We only expect 
the gas to have a velocity close to this estimate if it fell in 
freely from very large distances and if the Hubble expan- 
sion, which damps peculiar velocities, were unimportant. 
However, we do expect the scaling with mass to be more 
generally applicable. For the cold mode the trend with halo 
mass is indeed well reproduced by the escape velocity. 

On average, for halo masses above 10 11 Mg, the net 
accretion rate is positive, which means that more mass is 
flowing in than is flowing out (bottom- middle panel). There- 
fore, unsurprisingly, the haloes are growing. For haloes with 
10 10 M Q < M halo < 10 11 M Q the mean accretion rate is neg- 
ative (indicating net outflow), but small (~ 0.1 M^yr -1 ). 
The mean accretion rate of cold-mode gas is positive for 
all halo masses, but for hot-mode gas there is net outflow 
for Mhaio < 10 11 ' 5 Mq. Although these haloes are losing 
gas that is currently hot-mode, their hot-mode gas reser- 
voir may still be increasing if cold-mode gas is converted 
into hot-mode gas. For higher-mass haloes, the hot-mode ac- 
cretion rate is also positive and it increases approximately 
linearly with halo mass. This is the regime where the im- 
plemented supernova feedback is not strong enough to blow 
gas out of the halo. This transition mass is increased by 
more than an order of magnitude when more effective su- 
pernova feedback or AGN feedback is included (not shown). 
For Mhaio > 10 12 ' 5 M the hot-mode inflow rate is slightly 
stronger than the cold-mode inflow rate. 

The grey, dashed curve indicates the accretion rate a 
halo with a baryon fraction flb/n m would need to have to 
grow to its current baryonic mass in a time equal to the age 
of the Universe at z = 2, 



inverse 

Comparing this analytic estimate with the actual mean 
accretion rate, we see that they are equal for Mhaio > 
10 Mq, indicating that these haloes are in a regime of 
efficient growth. For lower-mass haloes, the infall rates are 
much lower, indicating that the growth of these haloes has 
halted, or that their baryon fractions are much smaller than 

fib / ftn ■ 

The bottom-right panel of Fig. [3] shows that, at the 
virial radius, hot-mode gas dominates the gas mass for high- 
mass haloes. The hot fraction at i? v ir increases from 10 per 
cent in haloes of ~ 10 10 M to 90 per cent for Mhaio ~ 
10 13 M©. Note that for haloes with Mhalo < 10 11,3 M the 
virial temperatures are lower than our adopted threshold 
for hot-mode gas. The hot fraction would have been much 
lower without supernova feedback (/hot < 5 per cent for 
Mhaio < lO 10 ' 5 M s ). 



5 INFLOW AND OUTFLOW 

Figs. [5] and [6] show the physical properties of the gas, 
weighted by the radial mass flux, for all, inflowing, and out- 
flowing gas (black, blue, and red curves, respectively). Ex- 
cept for the last two panels, the curves indicate the medians, 
i.e. half the mass flux is due to gas above the curves. Sim- 
ilarly, the shaded regions indicate the 16th and 84th per- 
centiles. Note that we do not plot this separately for hot- 
and cold-mode gas. The differences between the blue and 



red curves arise purely from the different radial peculiar ve- 
locity directions. 

We can immediately see that separating gas according 
to its radial velocity direction yields similar results as when 
the gas is separated according to its maximum past tem- 
perature (compare to Figs. [3] and [4]). Like cold-mode gas, 
inflowing gas has, on average, a higher density, a lower tem- 
perature, a lower entropy, and a lower metallicity than out- 
flowing gas. 

Comparing Figs. [3] and we notice that the differences 
in density, temperature, pressure, entropy, and metallicity 
between in- and outflowing gas tend to be slightly smaller 
than between cold- and hot-mode gas. This is particularly 
true outside the haloes (at R > 3i? v ir), where there is a clear 
upturn in the density and pressure of outflowing gas that 
is accompanied by a marked decrease in the temperature. 
These features are due to gas that is flowing towards other 
haloes and/or large-scale filaments. Although such gas is 
outflowing from the perspective of the selected halo, it is 
actually infalling gas and hence more likely to be cold-mode. 

Unsurprisingly, the radial peculiar velocities (bottom- 
left panels) are clearly very different when we separate the 
gas into in- and outflowing components than when we divide 
it into cold and hot modes. Low values of the radial velocity 
are avoided because the plot is weighted by the mass flux. 
While the radial velocity of cold-mode gas decreases from 
the virial radius towards the centre, the mass flux-weighted 
median radial velocity of inflowing gas is roughly constant. 
Within the haloes, the mass flux-weighted median radial ve- 
locity of outflowing gas decreases with radius. 

Both the inflow and the outflow mass flux (bottom- 
middle panel of Fig. [SJ are approximately constant inside 
0.7-Rvir, which implies that the fraction of the gas that is out- 
flowing is also constant (last panel of Fig. [5}. A mass flux 
that is independent of radius implies efficient mass trans- 
port, as the same amount of mass passes through each shell 
per unit of time. The inflowing mass flux decreases from 
10-Rvir to ~ i?vir because the overdensity of the region is 
increasing and because some of the gas that is infalling at 
distances S> -Rvir is falling towards neighbouring haloes. At 
R > 7? v ir the outflowing mass flux decreases somewhat, indi- 
cating that the transportation of outflowing material slows 
down and that the galactic winds are becoming less efficient. 
This can also be seen by the drop in outflow fraction around 
the virial radius in the last panel of Fig. [5] (The small de- 
crease in outflow fraction below 10 10 ' 5 Mq is a resolution 
effect.) The outflowing mass flux increases again at large 
radii, because the hot-mode gas is falling towards unrelated 
haloes. 

Comparing Figs. [4] and [6l we see again that, to first or- 
der, inflowing and outflowing gas behave similarly as cold- 
mode and hot-mode gas, respectively. There are, however, 
some clear differences, although we do need to keep in mind 
that some are due to the fact that Fig. [4] is mass- weighted 
while Fig. [S] is mass flux- weighted. The mass flux-weighted 
median temperature of outflowing gas is always close to the 
virial temperature. For Mhaio < 10 Mq this is much lower 
than the median temperature of hot-mode gas, but that 
merely reflects the fact that for these haloes the virial tem- 
perature is lower than the value of 10 5 ' 5 K that we use to 
separate the cold and hot modes. The mass flux-weighted 
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Figure 5. Properties of gas in haloes with 10 11 ' 5 Mq < Afhalo < 10 12 ' 5 Mq at z = 2 as a function of radius for all (black curves), 
outflowing (red curves), and inflowing (blue curves) gas. Results are shown for the simulation REF_L050N512. Shaded regions show 
values within the 16th and 84th percentiles, i.e., the ±l<x scatter around the median. From the top-left to the bottom-right, the different 
panels show the mass flux-weighted median gas overdensity, temperature, maximum past temperature, pressure, entropy, metallicity, 
radial peculiar velocity, the mean accretion rate, and the mean mass fraction of outflowing gas, respectively. Most of the trends with 
radius for inflowing and outflowing gas are similar to those for cold-mode and hot-mode gas, respectively, as shown in Fig. [3] 



maximum past temperature is about 0.5 dex higher than 

T- vir ■ 

Another clear difference is visible at the high-mass end 
(Mhaio > 10 12 ' 5 Mq). While the cold-mode density increases 
rapidly with mass, the overdensity of infalling gas remains ~ 
10 2 , and while the cold-mode temperature remains ~ 10 4 K, 
the temperature of infalling gas increases with halo mass. 
Both of these differences can be explained by noting that, 
around the virial radius, hot-mode gas accounts for a greater 
fraction of the infall in higher mass haloes (see the bottom- 
middle panel of Fig. 

As was the case for the cold-mode gas, the radial pe- 
culiar velocity of infalling gas scales like the escape veloc- 
ity (bottom- middle panel). Interestingly, although the mass 
flux- weighted median outflowing velocity is almost indepen- 
dent of halo mass, the high- velocity tail is much more promi- 
nent for low-mass haloes. Because the potential wells in these 
haloes are shallow and because the gas pressure is lower, the 
outflows are not slowed down as much before they reach the 
virial radius. The flux- weighted outflow velocities are larger 
than the inflow velocities for M^aio < 10 Mq, whereas 
the opposite is the case for higher-mass haloes. 

Finally, the last panel of Fig. [6] shows that the fraction 
of the gas that is outflowing around R v [ r is relatively stable 
at about 30-40 per cent. Although the accretion rate is neg- 
ative for 10 10 M Q < Mhaio < 10 11 M , which indicates net 
outflow, less than half of the gas is outflowing. 



6 EFFECT OF METAL-LINE COOLING AND 
OUTFLOWS DRIVEN BY SUPERNOVAE 
AND AGN 

Fig. [7] shows images of the same 10 12 Mq halo as Fig. [5] 
for five different high-resolution (L025N512) simulations at 
z = 2. Each row shows a different property, in the same order 
as the panels in the previous figures. Different columns show 
different simulations, with the strength of galactic winds in- 
creasing from left to right (although the winds are some- 
what stronger in NOZCOOL than in REF). The first col- 
umn shows the simulation without SN feedback and without 
metal-line cooling. The second column shows the simula- 
tion with SN feedback, but without metal-line cooling. The 
third column shows our reference simulation, which includes 
both SN feedback and metal-line cooling. The fourth column 
shows the simulation with density-dependent SN feedback, 
which is more effective at creating galactic winds for this 
halo mass. The last column shows the simulation that in- 
cludes both SN and AGN feedback. 

The images show substantial and systematic differences. 
More efficient feedback results in lower densities and higher 
temperatures of diffuse gas and in the case of AGN feedback, 
even some of the cold filaments are partially destroyed. Al- 
though the images reveal some striking differences, Fig. [8] 
shows that the trends in the profiles of the gas properties, 
including the differences between hot and cold modes, are 
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Figure 6. Properties of gas at 0.8-f? v ir < R < R v ii at z = 2 as a function of halo mass for all (black curves), outflowing (red curves), 
and inflowing (blue curves) gas. Results are shown for the simulation REF_L050N512. Shaded regions show values within the 16th and 
84th percentile, i.e., the ±lcr scatter around the median. From the top-left to the bottom-right, the different panels show the mass 
flux-weighted median gas overdensity, temperature, maximum past temperature, pressure, entropy, metallicity, radial peculiar velocity, 
the mean accretion rate, and the mean mass fraction of outflowing gas, respectively. The horizontal, dotted line in the first panel indicates 
the threshold for star formation (njj = 0.1 cm -3 ). The dashed, grey curves show analytic estimates from virial arguments. Most of the 
trends with mass for inflowing and outflowing gas are similar to those for cold-mode and hot-mode gas, respectively, as shown in Fig. [4] 



very similar in the different simulations. This is partly be- 
cause the profiles are mass-weighted, whereas the images 
are only mass-weighted along the projected dimension. We 
would have seen larger differences if we had shown volume- 
weighted profiles, because the low-density, high-temperature 
regions, which are most affected by the outflows, carry very 
little mass, but dominate the volume. Although the conclu- 
sions of the previous sections are to first order independent 
of the particular simulation that we use, there are some in- 
teresting and clear differences, which we shall discuss below. 

We can isolate the effect of turning off SN-driven out- 
flows by comparing models NOSN.NOZCOOL and NOZ- 
COOL. Without galactic winds, the cold-mode densities and 
pressures are about an order of magnitude higher within 
0.5i? v ir- On the other hand, turning off winds decreases the 
density of hot-mode gas by nearly the same factor around 
0.1i? v ir- Galactic winds thus limit the build up of cold-mode 
gas in the halo center, which they accomplish in part by 
converting cold-mode gas into hot-mode gas (i.e. by shock- 
heating cold-mode gas to temperatures above 10 5 ' 5 K). We 
can see that this must be the case by noting that the hot- 
mode accretion rate is negative inside the haloes for model 
NOZCOOL. The absence of SN-driven outflows also has a 
large effect on the distribution of metals. Without winds, the 
metallicities of both hot- and cold- mode gas outside the halo 
are much lower, because there is no mechanism to transport 



metals to large distances. On the other hand, the metallicity 
of the cold-mode gas is higher within the halo, because the 
star formation rates, and thus the rates of metal produc- 
tion, are higher. This suggests that the metals in cold-mode 
gas are associated with locally formed stars, e.g. infalling 
companion galaxies. 

Increasing the efficiency of the feedback, as in models 
WDENS and particularly AGN, the cold-mode median ra- 
dial velocity becomes less negative and the accretion rate of 
cold-mode gas inside haloes decreases. At the same time, the 
radial velocity and the outflow rate of the hot-mode increase. 
In other words, more efficient feedback hinders the inflow 
rate of cold-mode gas and boosts the outflows of hot-mode 
gas. The differences are particularly large outside the haloes. 
Whereas the moderate feedback implemented in model REF 
predicts net infall of hot-mode gas, the net accretion rate of 
hot-mode gas is negative out to about 4R v i r when AGN feed- 
back is included. Beyond the virial radius, stronger winds 
substantially increase the mass fraction of hot-mode gas. 
This comes at the expense of the hot-mode gas inside the 
haloes, which decreases if the feedback is more efficient, at 
least for 0.2i? v ir < R < Rviv 

The effect of metal-line cooling can be isolated by com- 
paring models NOZCOOL and REF. Without metal-line 
cooling, the cooling times are much longer. Consequently, 
the median temperature of the hot-mode gas remains above 
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Figure 7. As Fig. [2] but comparing the different simulations listed in Table [2] The images show a cubic 1 h^ 1 comoving Mpc region 
centred on a halo of M halo RJ 10 12 Mq at z = 2 taken from the L025N512 simulations. From left to right, the different columns show the 
simulation with neither SN feedback nor metal-line cooling, the simulation without metal-line cooling, the reference simulation (which 
includes SN feedback and metal-line cooling), the simulation for which the wind speed scales with the effective sound speed of the ISM, 
and the simulation with AGN feedback. Note that the images shown in the middle column are identical to those shown in Fig. [2] 
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Figure 8. Properties of gas in haloes with 10 11 ' 5 Mq < Af^alo < 10 12,5 Mq at z = 2 as a function of radius for hot-mode (red curves) 
and cold-mode (blue curves) gas for the five different L025N112 simulations listed in Table |21 as indicated in the legend. From the top- 
left to the bottom-right, the different panels show the mass-weighted median gas overdensity, temperature, maximum past temperature, 
pressure, entropy, metallicity, radial peculiar velocity, the mean accretion rate, and the mean mass fraction of hot-mode gas, respectively. 



10 6 K (at least for R > 0.17? v ir), whereas it suddenly drops 
to below 10 K around 0.2i? v ir when metal-line cooling is 
included. Thus, the catastrophic cooling flow of the diffuse, 
hot component in the inner haloes is due to metals. Indeed, 
while the median hot-mode radial peculiar velocity within 
0.27? v i r is positive without metal-line cooling, it becomes 
negative (i.e. infalling) when metal-line cooling is included. 



7 EVOLUTION: MILKY WAY-SIZED HALOES 

AT Z = 

Fig. [5] is identical to Fig. [3] except that it shows profiles 
for z = rather than z — 2. For comparison, the dotted 
curves in Fig. [9] show the corresponding z = 2 results. As 
we are again focusing on 10 115 M < M halo < 10 125 Mq, 
the results are directly relevant for the Milky Way galaxy. 

Comparing Fig. [3] to Fig. [9] (or comparing solid and 
dotted curves in Fig. 0), we see that the picture for z — 
looks much the same as it did for z = 2. There are, however, 
a few notable differences. The overdensity profiles hardly 
evolve, although the difference between the cold and hot 
modes is slightly smaller at lower redshift. However, a con- 
stant overdensity implies a strongly evolving proper density 
(p oc (1 + z) 3 ) and thus also a strongly evolving cooling rate. 
The large decrease in the proper density caused by the ex- 
pansion of the Universe also results in a large drop in the 
pressure and a large increase in the entropy. 

The lower cooling rate shifts the peak of the cold-mode 



temperature profile from about 2i? v ir to about 7? v ir- While 
there is only a small drop in the temperature of the hot- 
mode gas, consistent with the mild evolution of the virial 
temperature of a halo of fixed mass (T V i r oc (1 + z); eq. [7]), 
the evolution in the median temperature for all gas is much 
stronger than for the individual accretion modes. While at 
z = 2 the overall median temperature only tracks the median 
hot-mode temperature around 0.5_R v ir < R < 2J? v ir, at z — 
the two profiles are similar at all radii. 

The metallicity profiles do not evolve much, except for 
a strong increase in the median metallicity of cold-mode gas 
at R S> Rvir- Both the weak evolution of the metallicity of 
dense gas and the stronger evolution of the metallicity of the 
cold, low-density intergalactic g as far away from galax ies are 
consistent with the findings of W lersma et alJ (|201ll ). who 
found these trends to be robust to changes in the subgrid 
physics. 

The absolute, net radial velocities of both the hot- 
and cold-mode components are smaller at lower redshift, 
as expected from the scaling of the characteristic velocity 
(v csc oc (l + z) 1/2 ; eq. QT]). 

While at z — 2 the net accretion rate was higher for 
the cold mode at all radii, at z — the hot mode domi- 
nates beyond 37? v i r . At low redshift the rates are about an 
order of magnitude lower than at z = 2. For R < 0.4i? v i r 
the net accretion rate is of order 1 Mq yr _1 , which is domi- 
nated by the cold mode, even though most of the mass is in 
the hot mode. Since a substantial fraction of both the cold- 
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Figure 9. Properties of gas in haloes with 10 11 ' 5 Mq < M^aio < 10 12 ' 5 Mq at z = as a function of radius for all (black curves), 
hot-mode (red curves), and cold-mode (blue curves) gas. Results are shown for the simulation REF_L050N512. Most z = 2 curves from 
Fig. [3] have also been included, for comparison, as dotted curves. Shaded regions show values within the 16th and 84th percentiles, i.e., 
the ±1<t scatter around the median. From the top-left to the bottom-right, the different panels show the mass-weighted median gas 
overdensity, temperature, maximum past temperature, pressure, entropy, metallicity, radial peculiar velocity, the mean accretion rate, 
and the mean mass fraction of hot-mode gas, respectively. The radial profiles at z = follow the same trends as at z = 2 (compare 
Fig. [3} , although the pressure, the infall velocity and the accretion rates are lower, the entropy is higher, and the hot mode accounts for 
a greater fraction of the mass. 



and hot-mode gas inside this radius is outflowing, the actual 
accretion rates will be a bit higher. 

At z — the fraction of the mass that has been hotter 
than 10 5 ' 5 K exceeds 50 per cent at all radii and the profile 
show a broad peak of around 80 per cent at 0.3 < R < R v i T . 
Thus, hot-mode gas is more important at z = than at 
z = 2, where cold-mode gas accounts for most of the mass 
for R < 0.3R v j T and R > 27? v i r . This is consistent with 
Ivan de Voort et al.l (|2011al ). who investigated the evolution 
of the hot fraction in more detail using the same simulations. 



8 CONCLUSIONS AND DISCUSSION 

We have used cosmological hydrodynamical simulations 
from the OWLS project to investigate the physical prop- 
erties of gas in and around haloes. We paid particular at- 
tention to the differences in the properties of gas accreted 
in the cold and hot modes, where we classified gas that has 
remained colder (has been hotter) than 10 K while it was 
extragalactic as cold-mode (hot-mode) gas. Note that our 
definition allows hot-mode gas to be cold, but that cold- 
mode gas cannot be hotter than 10 5 ' 5 K. We focused on 
haloes of 10 12 M at z = 2 drawn from the OWLS ref- 
erence model, which includes radiative cooling (also from 



heavy elements), star formation, and galactic winds driven 
by SNe. However, we also investigated how the properties of 
gas near the virial radius change with halo mass, we com- 
pared z = 2 to 2 = 0, we measured properties separately 
for inflowing and outflowing gas, and we studied the effects 
of metal-line cooling and feedback from star formation and 
AGN. We focused on mass-weighted median gas properties, 
but noted that volume-weighted properties are similar to 
the mass-weighted properties of the hot-mode gas, because 
most of the volume is filled by dilute, hot gas, at least for 
haloes with T vir > 10 5 5 K (see Figs, [l] and [2]). 

Let us first consider the properties of gas just inside 
the virial radius of haloes drawn from the reference model 
at z = 2 (see Fig. |4}. The fraction of the gas accreted in 
the hot mode increases from 10 per cent for halo masses 
Mhalo ~ 10 10 M Q to 90 per cent for M ha io ~ 10 13 M . 
Hence, 10 12 M is a particularly interesting mass scale, as 
it marks the transition between systems dominated by cold 
and hot mode gas. 

Although the cold streams are in local pressure equilib- 
rium with the surrounding hot gas, cold-mode gas is physi- 
cally distinct from gas accreted in the hot mode. It is colder 
(T < 10 5 K vs. T > T vir ) and denser, particularly for high- 
mass haloes. While hot-mode gas at the virial radius has a 
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density ~ 10 2 (p) for all halo masses, the median density of 
cold-mode gas increases steeply with the halo mass. 

While the radial peculiar velocity of cold-mode gas is 
negative (indicating infall) and scales with halo mass like the 
escape velocity, the median hot-mode velocities are positive 
(i.e. outflowing) for Mhaio <■ 10 M© and for larger masses 
they are much less negative than for cold-mode gas. Except 
for Mhaio ~ 10 10 ' 5 Mq, the net accretion rate is positive. For 
10 12 < A/haio < 10 13 Mq the cold- and hot-mode accretion 
rates are comparable. 

While hot-mode gas has a metallicity ~ 10 _1 Zq, the 
metallicity of cold-mode gas is typically significantly smaller 
and displays a much larger spread. The scatter in the local 
metallicity of cold-mode gas is large because the cold fila- 
ments contain low-mass galaxies that have enriched some of 
the surrounding cold gas. We emphasized that we may have 
overestimated the median metallicity of cold-mode gas be- 
cause we find it to decrease with increasing resolution. It is 
therefore quite possible that most of the cold-mode gas in 
the outer halo still has a primordial composition. 

The radial profiles of the gas properties of haloes with 
10 11 ' 5 < A/ h aio < lO 12 ' 5 M Q revealed that the differences 
between gas accreted in the cold and hot modes vanishes 
around 0.1i? v ir (see Fig. [3]), although the radius at which this 
happens decreases slightly with the resolution. The conver- 
gence of the properties of the two modes at small radii is due 
to catastrophic cooling of the hot gas at R < 0.2i? v i r . Inter- 
estingly, in the absence of metal-line cooling the hot-mode 
gas remains hot down to much smaller radii, which suggests 
that it is very important to model the small-scale chemical 
enrichment of the circumgalactic gas. While stronger winds 
do move large amounts of hot-mode gas beyond the virial ra- 
dius, even AGN feedback is unable to prevent the dramatic 
drop in the temperature of hot-mode gas inside 0.2_R v ir, at 
least at z = 2. 

While the density and pressure decrease steeply with ra- 
dius, the mass-weighted median temperature peaks around 
0.5 — 1.0-Rvir- This peak is, however, not due to a change in 
the temperature of either the cold- or hot-mode gas, but due 
to the radial dependence of the mass fraction of gas accreted 
in the hot mode. The hot-mode fraction increases towards 
the halo, then peaks around 0.5 — 1.0i? v ir and decreases mov- 
ing further towards the centre. Even outside the halo there 
is a significant amount of hot-mode gas (e.g. ~ 30 per cent 

at 10i?vir). 

Beyond the cooling radius, i.e. the radius where the 
cooling time equals the Hubble time, the temperature of the 
hot-mode gas decreases slowly outwards, but the tempera- 
ture of cold-mode gas only peaks around 2R v i r at values just 
below 10 5 K. The density for which this peak temperature is 
reached depends on the interplay between cooling (both adi- 
abatic and radiative) and photo-heating. The metallicity de- 
creases with radius and does so more strongly for cold-mode 
gas. Near 0.1i? v ir the scatter in the metallicity of cold-mode 
gas is much reduced, although this could be partly a resolu- 
tion effect. 

The median radial peculiar velocity of cold-mode gas is 
most negative around 0.5 — 1.0i? v ir- For the hot mode, on 
the other hand, it is close to zero around that same radius. 
Hence, the infall velocity of the cold streams peaks where the 
hot-mode gas is nearly static, or outflowing if the feedback 
is very efficient. We note, however, that the scatter in the 



peculiar velocities is large. For R < J? v j r much of the hot- 
mode gas is outflowing and the same is true for cold-mode 
gas at R ~ 0.1-Rvir- 

Inside the halo the cold-mode accretion rate increases 
only slightly with radius (d In M/d In R ~ 0.4), indicating 
that most of the mass is transported to the central galaxy. 
For the hot-mode, on the other hand, the accretion rate only 
flattens at R < 0.4i? v ir- This implies that the hot accretion 
mode mostly feeds the hot halo. However, hot-mode gas that 
reaches radii ~ 0.17? v ir is efficiently transported to the centre 
as a result of the strong cooling flow. Nevertheless, cold- 
mode accretion dominates the accretion rate at all radii. 

Dividing the gas into inflowing and outflowing compo- 
nents yields results that are very similar to classifying the 
gas on the basis of its maximum past temperature. This is 
because inflowing gas is mostly cold-mode and outflowing 
gas is mostly hot-mode. The situation is, however, different 
for high- mass haloes (Mhaio > 10 12 ' 5 Mq). Because the two 
accretion modes bring similar amounts of mass into these 
haloes, the properties of the infalling gas are intermediate 
between those of the cold and hot accretion modes. 

When expressed in units of the mean density of the 
universe, the z — density profiles (with radius expressed 
in units of the virial radius) are very similar to the ones 
at z = 2. The same is true for the metallicities and tem- 
peratures, although the peak temperatures shift to slightly 
smaller radii (again normalized to the virial radius) with 
time. A fixed overdensity does imply that the proper den- 
sity evolves as p oc (1 + z) 3 , so the pressure (entropy) are 
much lower (higher) at low redshift. Infall velocities and ac- 
cretion rates are also significantly lower, while the fraction 
of gas accreted in the hot mode is higher. The difference in 
the behaviour of the two accretion modes is, however, very 
similar at z — and z = 2. 

Although there are some important differences, the 
overall properties of the two gas modes are very similar be- 
tween our different simulations, and are therefore insensitive 
to the inclusion of metal-line cooling and galactic winds. 
Without SN feedback, the already dense cold-mode gas 
reaches even higher densities, the metallicity of hot-mode 
gas in the outer halo is much lower. Without metal-line cool- 
ing, the temperature at radii smaller than 0.2i? v ir is much 
higher. With strong SN feedback or AGN feedback, a larger 
fraction of the gas is outflowing, the outflows are faster, and 
the peak in the fraction of the gas that is accreted in the hot 
mode peaks at a much larger radius (about 3-R v ir instead of 
0.5 - l.Oflvh-). 

iKeres et all l|201ll) have recently shown that some of 
the hot gas properties depend on the numerical technique 
used to solve the hydrodynamics. They did, however, not 
include metal-line cooling and feedback. They find that the 
temperatures of hot gas (i.e. T > 10 5 K) around 10 12-13 M Q 
haloes are the same between the two methods, but that the 
median density and entropy of 10 12 Mq haloes is somewhat 
different, by less than a factor of two. These differences are 
comparable to the ones we find when using different feedback 
models. They also show that the radial velocities are differ- 
ent at small radii. Our results show that including metal-line 
cooling decreases the radial velocities, whereas including su- 
pernova or AGN feedback increases the radial velocities sig- 
nificantly. Even though there are some differences, the main 
conclusions of this work are unchanged. The uncertainties 
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associated with the subgrid implementation of feedback and 
the numerical method are therefore unlikely to be important 
for our main conclusions. 

Cold (i.e. T <C 10 5 K) outflows are routinely detected 
in the form of blueshifted interstellar absorption lines in 
the rest-frame UV spectra of star-forming galaxies (e.g . 
Weiner et all 120091 ; Isteidel et all l20ld ; iRubin et all I201C 



i.e. A^hi > 10 172 cm 2 ) may b e used to trace co l d flows 



Rakic et al.ll2011af i. This is not in conflict with our results, 



because we found that the outflowing gas spans a very wide 
range of temperatures and because the detectable UV ab- 
sorption lines are biased towards colder, denser gas. Addi- 
tionally, the results we showed are mass- weighted, but if we 
had shown volume- weighted quantities, the outflow fractions 
would have been larger. The inflowing material has smaller 
cross-sections and is therefore less likely to be detec ted (e.g. 
IFaucher-Giguere fc Keresj|201ll ; IStewart et al.ll201l| y 

How can we identify cold-mode accretion observation- 
ally? The two modes only differ clearly in haloes with 
T vir > f0 5 K (M halo » fO 10 - 5 M ((l + 2)/3)~ 3/2 ), because 
photo-ionization by the UV background radiation ensures 
that all accreted gas is heated to temperatures up to ~ f 5 K 
near the virial radius. Near the central galaxy, R < 0.1-R v ir, 
it is also difficult to distinguish the two modes, because gas 
accreted in the hot mode is able to cool. 

In the outer parts of sufficiently massive haloes the 
properties of the gas accreted in the two modes do differ 
strongly. The cold-mode gas is confined to clumpy filaments 
that are approximately in pressure equilibrium with the dif- 
fuse, hot-mode gas. Besides being colder and denser, cold- 
mode gas typically has a much lower metallicity and is much 
more likely to be infalling. However, the spread in the prop- 
erties of the gas is large, even for a given mode and a fixed 
radius and halo mass, which makes it impossible to make 
strong statements about individual gas clouds. Neverthe- 
less, it is clear that most of the dense (p S> 10 2 (p)) gas in 
high-mass haloes (T v i r > fO 6 K) is infalling, has a very low 
metallicity and was accreted in the cold-mode. 

Cold-mode gas could be observed in UV line emis- 
sion if we are able to detect it in the outer halo of 
massive galaxies. Diffuse Lyman-a emission has already 
been detected (e.g. ??), but its interpretation is compli- 
cated by radiative transfer effects and the detected emis- 
sion is mor e likely scattere d ligh t from central Hll re- 
gions (e .g. iFurlanetto et alJ 120051: IFaucher-Giguere et all 



20ld: ?; lHaves et al. 



2011 



iRauch et alJ 12011 



_ _ a, but see also 

Diikstra fc Loebll200gl : iRosdahl fc Blaizodl201lh . Metal- line 
emission from ions such as C III, C IV, Si III, and Si IV could 
potentially reveal cold streams, but current facilities do not 
probe down to the expected surface brightnesses (??). 

The typical temperatures (T ~ 10 4 K) and densities 
(P ^ 10 2 (p)) correspond to those of strong quasar absorp- 
tion lines systems. For example, at z — 2 the typ ical Hi 
colum n density is0 JV H , ~ 10 16 cm~ 2 (p/[10 2 (p>]) 3/2 l|Schavel 
l2001bl ) with higher column density gas more likely to have 
been accreted in the cold mode. At low redshift the H I col- 
umn densities corresponding to a fixed overden sity are about 
1-2 orders of magnitude lower (|Schavdl2001bl ). 

Indeed, simulations show that Lyman limit systems 



3 For Nwi > 10 18 cm " 2 
see e.g. ISchavdEuOlal : lAltav et al]|201ll . 



Faucher-Giguere fc Kerel l201ll: iFumagalli et alJ 12011; 
Ivan de Voort et alJl2oiTd ) and Ivan de Voort et all l|2011(J) 
have demonstrated that cold-mode accretion is required to 
match the observed rate of incidence of strong absorbers 
at z = 3. Many stro ng QSO absorbers also tend to have 
low metallicities (e.g. iRibaudo et alJl201ll : iGiavalisco et al.l 
l201ll : IFumagalli et~aj]l2011bh . although it should be noted 
that metallicity measurements along one dimension may un- 
derestimate the mean metallicities of three-dimensional gas 
clouds due to the e xpected poor small-scale metal mixing 
jSchave et all 120071 ). We also note that most Lyman limit 
systems are predicted to arise in or around haloes with 
masses that are much lower than required for the presence of 
stable accretion shocks near the virial radius, so that they 
will generally not corr espond to cold streams p enetrating 
hot, hydrostatic haloes (|van de Voort et al.ll2011q |. To study 
those, it is therefore more efficient to target s ight lines to 



QSOs close to massive foreground galaxies (e.g. Rakic etaD 



2011bl: IFaucher-Giguere fc Keresl I20TTI : IStewart et aDl201ll : 



Kimm et aUl201lr ). 
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APPENDIX A: RESOLUTION TESTS 

We have checked (but do not show) that the results pre- 
sented in this work are converged with respect to the size of 
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Figure Al. Convergence of gas profiles with resolution. All profiles are for haloes with 10 11 ' 5 M < M halo < 10 12 - 5 M at z = 2. 
Results are shown for simulations REF.L025N512 (solid), REF.L050N512 (dotted), and REF.L050N512 (dashed) and for both the 
hot (red) and cold (blue) modes. The three simulations vary by a factor of 64 (8) in mass (spatial) resolution. From the top-left to the 
bottom-right, the different panels show the mass-weighted median gas overdensity, temperature, maximum past temperature, pressure, 
entropy, mctallicity, radial peculiar velocity, the mean accretion rate, and the mean mass fraction of hot-mode gas, respectively. 




Figure A2. Convergence of the properties of gas at 0.8-R v i r < R < fi? v ir as a function of halo mass with resolution. Panels and curves 
as in Fig. lAf I 
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the simulation volume if we keep the resolution fixed. The 
only exception is hot-mode gas at R > 5R v i T , for which the 
radial velocities and accretion rate require a box of at least 
50/i _1 Mpc on a side (which implies that our fiducial simu- 
lation is sufficiently large). 

Convergence with resolution is, however, more difficult 
to achieve. In Figs. I All and IA2I we show again the radial 
profiles and mass dependence of the halo properties for the 
hot- and cold-mode components at z — 2. Shown are three 
different simulations of the reference model, which vary by a 
factor of 64 (8) in mass (spatial) resolution. All trends with 
radius and halo mass are very similar in all runs, proving 
that most of our conclusions are robust to changes in the 
resolution. Below we will discuss the convergence of hot- 
and cold-mode gas separately and in more detail. 

The convergence is generally excellent for hot mode gas. 
As the resolution is increased, the density of hot-mode gas 
decreases slightly and the temperature drop close to the halo 
centre shifts to slightly smaller radii, which also affects the 
pressure and entropy. There is a small upturn in the den- 
sity of hot-mode gas at the virial radius as we approach the 
halo mass corresponding to the imposed minimum of 100 
dark matter particles, showing that we may have to choose 
a minimum halo mass that is a factor of 5 higher for complete 
convergence. The radial peculiar velocity increases slightly 
with resolution, causing the hot-mode accretion rate to de- 
crease. 

Convergence is more difficult to achieve for cold-mode 
gas. The density of cold-mode gas inside haloes, and thereby 
also the difference between the two modes, increases with 
the resolution. The pressure of the cold-mode gas also in- 
creases somewhat with resolution, which leads to a smaller 
difference with the pressure of hot-mode gas. The cold-mode 
radial velocity becomes more negative, increasing the differ- 
ence between the two modes. 

The median metallicity of the cold-mode gas decreases 
strongly with increasing resolution for R > 0.3i? v ir- The 
metallicity difference between the two modes therefore in- 
creases, although the distributions still overlap (not shown), 
and the radius at which the metallicities of the two modes 
converge decreases. In fact, the convergence of the median 
metallicity of cold-mode gas is so poor that we cannot rule 
out that it would tend to zero at all radii if we keep increas- 
ing the resolution. 

If we use particle metallicities rather than SPH 
smoothed metallicities (see Section |3.6|) . then the median 
metallicity of both modes is lower and the median cold-mode 
metallicity plummets to zero around the virial radius (not 
shown). The decrease in metallicity with increasing resolu- 
tion is in that case less strong, but still present. The un- 
smoothed hot-mode metallicities are also not converged and 
lower than the (converged) smoothed hot-mode metallicities, 
but they increase with increasing resolution. The difference 
between the two modes therefore increases with resolution. 

With increasing resolution, the radial velocities of cold- 
mode gas become slightly more negative within the halo. 
The net accretion rates and hot fractions are converged. 

Even though some properties are slightly resolution de- 
pendent, or strongly so for the case of the metallicity of 
cold-mode gas, all our conclusions are robust to increases in 
the numerical resolution. 
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